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Abstract 

The transport phenomena of a nonequiUbrium lattice gas system are investigated. We consider a 
simple system that consists of two particles interacting repulsively and the potential forces acting 
on these particles. Under an external driving field applied to only one particle, we found the 
following relation between the mean velocity of the driven particle and the coefficient of effective 
drag of this particle under certain conditions; With the increase in the mean velocity, the coefficient 
of effective drag varies in the form, increase — > decrease — > increase — > decrease. Moreover, under 
other conditions, we found the following relations between these values which show changes similar 
to those between shear rate and shear viscosity observed in the shear-thickening polymer or colloidal 
suspensions; With the increase in the mean velocity, the coefficient of effective drag varies in the 
form, increase — > decrease or decrease — > increase — > decrease. We explain the mechanisms of such 
phenomena by considering the transition diagrams. 
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INTRODUCTION 



Nonequilibrium lattice gases are simple mathematical models, which have been useful 
and important in studies of the several properties of nonequilibrium systems with numer- 
ous degrees of freedom^]. Recently, many varieties of nonequilibrium phenomena such as 

appearances of 



nonequilibrium phase transitions 2 1 , a variety of particle flows s], 0, 0, Q , 

8|, and the mathematical foundations of nonequilibrium 



11 



12| | have been investigated through such 



long-range spatial correlations (7|, 
statistical mechanics and thermodynamics ol-liol 
lattice gases. 

In this paper, we investigate the transport behaviors of a lattice gas system with a periodic 
boundary, which consists of only two particles interacting repulsively and the potential forces 
acting on them. Under nonequilibrium conditions, lattice gas systems have been known to 
show some nontrivial phenomena, such as the appearance of long-range spatial correlations 0] 
and anomalous drift motions^], even if the system involves only two particles. 

In our system, the following complex transport properties are found when only one parti- 
cle is driven by an external driving field; With the increase in the mean velocity of the driven 
particle, the coefficient of effective drag of this particle (= [driving field strength] /[mean ve- 
locity]) varies in the form, increase — decrease increase — >■ decrease under certain con- 
ditions. Moreover, under other conditions, the relation between these values show changes 
similar to those between shear rate and shear viscosity observed experimentally in the shear- 
thickening polvmer jl^. Ilfil. Ilfi| or colloidal suspensions [21 . 22 1; With the increase in the mean 
velocity, the coefficient of effective drag varies in the form, increase — >■ decrease or decrease 

increase decrease. 

In the following sections, details of the numerical and analytical results are shown. First, 
we introduce our model and show the numerical results, where the above-mentioned phe- 
nomena are observed. Next, we explain the mechanism of these phenomena by considering 
the transition diagrams. Last, we show the similarities between our obtained results and 
those observed in the polymer solutions or colloidal suspensions. 
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MODEL 



Now, we introduce a lattice gas model, which is the same as that studied in our previous 
paper p. We consider a lattice system with two parallel one-dimensional lanes where each 
lane involves L sites with a periodic boundary. Each lane contains only one particle which 
moves randomly to the nearest sites without changing lanes. The sites occupied by particles 
in the 1st and 2nd lanes are denoted Xi and X2i respectively, which are given as integer 
numbers from to L — 1. 

The effect of potential forces acting on the particles is described by the following Hamil- 
tonian: 

H{xi, X2) = V{X,) + V{X2) + Vu{Xi,X2), (1) 

where V{x) represents the one-body potential on each lane, and V^i2(a^i, 2:2) represents the 
interaction potential between the two particles. Furthermore, an external driving field is 
applied to the particle on the 2nd lane. We denote the field strength F. 

The time evolution of this system is described by the iteration of the following three 
steps. First, one of the two particles is randomly chosen. Let the position of the chosen 
particle be x. Second, its neighboring site y, x — 1 or x + 1, is randomly chosen. Third, the 
chosen particle moves from x to y with the following probability 

cix,y, Xi,X2) I Qxp[Q{x ^ y;xi,X2)/kBTy ^^-^ 

with 

Q{x y; xi, X2) = H{x[, x'2) - H{xi, X2) - F{x2 - X2), (3) 

where {x[,X2) = {xi,y) when x = X2, and {x[,X2) = {y,X2) when x = xjj^. T is tempera- 
ture and the Boltzmann constant ks is set 1. Here, the time step is given by [No. of above 
iterations] /[No. of particles (= 2)]. 

Specifically, we study the case where V{x) = V\L/2—x\ (Fig. 1), and Vi2(a;i, X2) = ISxj^,x2 
using the L x L unit matrix 6ij. Also, we focus on the case L = 4. We found that this size is 
the minimum required to exhibit the phenomenon we demonstrate in the presented paper. 
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SIMULATION RESULTS 



In this section, we show the simulation results of this system. We mainly focus on the 
cases with V < I and \F\ < I+V because we obtained some novel phenomena as shown in the 
belows. Such conditions indicate that the influences of the potential forces and the driving 
field are weak compared to those of the particle-particle interactions. In order to characterize 
the system, we define the mean velocity of the driven particle (in the 2nd lane) in steady 
state, u, as the difference of the long time average of the moving ratio of the driven particle 
in the positive and negative directions. Here, the direction Xi :0— > 1— > > {L—1) — >0— > 
is positive. For simphcity, I — 1 and F > are set. 

First, we focus on the cases where T is small enough compared to /. Figures 2(a), (b) 
and (c) show m as a function of F for (a) V = 0.25 with T = 0.05 or T = 0.07, (b) V = 0.4 
with T = 0.05 or T = 0.08, and (c) V = 0.6 with T = 0.08 or T = 0.1. Then, the following 
two types of F — u profiles appear; A) u increases slowly —>■ steeply slowly steeply 
with the increase in F as shown in (a) and (b) with a smaller T. B) u increases steeply — > 
slowly — > steeply with the increase in F as shown in (b) with a larger T and (c). 

From these results, the relations between u and the coefficient of effective drag of the 
driven particle, r], defined as F/u are straightforwardly obtained. Figures 3(a), (b) and (c) 
show 1] as a function of u for (a) V = 0.25 ot V = 0.3, (b) V = 0.35 or V = 0.4, and (c) 
V = 0.55 or V = 0.6. Here, the properties oiu — r] profile change depending on V and T as 
following. 

I) -i) In the case with V < 7/3, as shown in Fig. 3(a), rj varies in the form, increase — > 
decrease increase — > decrease, with the increase in u for small T (for example T — 0.05). 
-ii) On the other hand for a little larger T (for example T = 0.07), the change in r] becomes 

less sharp and changes in the form, increase fast decrease slow decrease — fast decrease, 
with the increase in u. Here, u at the maximum values of rj decrease with the increase in V 
or the decrease in T. 

II) -i) In the case with 7/3 < V < 1/2, as shown in Fig. 3(b), rj varies in the form, 
increase — > decrease — > increase — > decrease, with the increase in u for small T (for example 
T — 0.05). -ii) On the other hand for a little larger T (for example T = 0.08), the change 

in 7] becomes less sharp, and simpler in the form, decrease increase decrease, with the 
increase in u. Here, u at the maximum values of r] decrease with the increase in V or the 
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decrease in T. 

Ill) In the case with 1/2 < V, as shown in Fig. 3(c), rj varies in the form, decrease — *• 
increase decrease, with the increase in u independently of T. Here, u at the maximum 
values of t] decrease with the decrease in V or T. 

It is remarkable that the u — 7] profile given in case I-ii) appears qualitatively similar to 
that between the shear rate and shear viscosity coefficient of the shear-thickening polymer 



solutions obtained experimentally 
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mm- 

On the other hand, in the case Il-ii) or III), the 



u — 7] profile appear qualitatively similar to that between the shear rate and shear viscosity 
coefficient of the shear-thickening colloidal suspension obtained experimentally^, 13]. 

Phase diagram for the properties oi u — 7] profile as functions of V and T is given as 
Fig. 3(d). Here, the boundaries given as = J/3 and V < I /2 (as indicated by the broken 
lines.) are predicted by the analytical studies as mentioned in the next section. If T is given 
larger, u — rj profiles for several V become less sharp and close to fiat like Newtonian fiuid 
continuously as shown in Fig. 3(a), (b) and (c) (for example T = 0.2 or T = 0.32). 

It is noted that, if \^ >> / in which the particle-particle interactions are much smaller 
than any other effects, u is roughly estimated as oc exp(— (2V" — / — F)/T) — exp(— (2V — 
I + F) /T). In such cases, 1] decreases monotonically with the increase in u for small T, and 
7] — u profile closes to fiat with the increase in T. Thus, the particle-particle interactions 
play important roles to show the presented non-trivial behaviors like those observed shear- 
thickening solutions. 



ANALYSIS OF THE MODEL BY TRANSITION DIAGRAMS 

In this section, we try to explain the mechanism of our obtained phenomena by considering 
the transition diagrams.. First, we name all states of this system using the sites occupied 
by the particles in each lane (xi,X2). Since L = 4, this system has 16 states from (0,0) 
to (3,3). The transitional tendency from {xi,X2) to is given by Q{x y;xi,X2), 

which was defined previously. Here, ( y) when x = X2, and {x[,X2) = (y, ^2) 

when X = Xi- Then, we obtain the transition diagrams; including all information on the 
transitions between states. 

We focus on the F dependent changes of the transition diagrams. Figure 4 shows the 
typical transition diagrams in the cases where (a) F < V, (h) V < F < I — 2V, (c) 
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I -2V < F < I -V, (d) I -V < F < I and (e) I < F under V < 1/3. Here, indicates 
the state (xi, X2), the arrows give the direction of the transition with a probabihty higher 
than ~ 0.5 implying Q{x y;xi,X2) ~ 0, and the values beside these arrows indicate 
the value of \Q{x y;xi,X2)\. The motion of the system is represented by the random 
walk caused by the fluctuation with T in the transition diagram. Here, the transition in a 
vertical direction indicates the motion of the particle in the 1st lane, and that in a horizontal 
direction indicates the motion of the driven particle in the 2nd lane. 

As shown in Figs. 4, there exist two types of states, the transient state, which has some 
arrows directing to some other states, and the attractor, which has no arrows directing to 
any other states. If T is small enough compared to / and V, the transitions from each 
state arc considered to conform to the following tendencies in most cases. 1) From the 
transient states, one of all transitions with the arrows directing to any other of the states 
is realized randomly. 2) From the attractors, only a transition in the direction with the 
smallest \Q{x — > y;xi,X2)\ {— \Qmin\) is realized with the probability ~ exp{—Qmin/T) 
(escape rate). It is noted that the system tends to stay the attractors for most of the time. 
Thus, the escape rate from these attractors contributes dominantly to give the mean velocity 
of the driven particle u. 

Now, we focus on the typical trajectories in the transition diagrams for the cases with 
V < 1/3 (Figs. 4) to explain the appearance of the F — u relations hke in Fig. 2(a). When 
F < V, the transition diagram as given in Fig. 4(a) involves the following trajectory as 
a typical one characterizing the steady state motion, (1,2) — > (1,3) — > (2,3) — > (2,0) — > 
(2,1) (3,1) (3,2) (0,2) (1,2). (as indicated by thin arrows). This trajectory 
indicates the motions in which both particles proceed in the direction of F . Here, this 
involves four attractors, (1,2), (2,3), (2, 1) and (3, 2). The escape rates from the attractors 
(1, 2) or (2, 3) given through the transitions (1,2) (1, 3) or (2, 3) (2,0) increase with the 
increase in F. However, those of (2,1) or (3, 2) given through the transitions (2, 1) — > (3, 1) 
or (3, 2) — > (0, 2) are independent of F. Then, u cannot increase so fast with F in this case. 

In the case where V < F < I — 2V, the typical trajectories giving the steady state 
motions still involve the attractors (2, 1) and (3,2) as shown in Fig. 4(b) (as indicated by 
thin arrows). Also in this case, the escape rates from these attractors are still independent 
of F. In this case, u increases little with F. 

In the case where I — 2V < F < I — V, the following two types of typical trajectories 
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appear, as shown in Fig. 4(c), i) trajectories including two attractors, (2, 1) and (3,2), for 
example (1,2) ^ (1,3) ^ (2,3) ^ (2,0) ^ (2,1) ^ (2,2) ^ (3,2) ^ (0,2) ^ (1,2) 
(as indicated by thin arrows), and ii) trajectories including only one attractor (2,1) for 
example (2,1) (2,2) (2,3) (2,0) — * (2,1) (as indicated by broken arrows). Here, 
the former trajectory indicates the motions in which both particles proceed in the direction 
of F while the latter indicates the motions in which only the driven particle proceeds. In 
this case, the escape rate from (2, 1) given through the transition (2, 1) (2, 2) increases 
exponentially with the increase in F while the escape rate (3, 2) given through the transition 
(3, 2) (0, 2) is independent of F. Then, the mean velocity of the driven particle remains 
unaffected by the increase in F if the system develops along trajectories i). On the other 
hand, this increases exponentially with F if the system develops along trajectories ii). Thus, 
within this range oi F, u increases steeply with F on the average. 

In the case where I — V < F < I, similar to the previous case, the following two types of 
trajectories are still typical ones as shown in Fig. 4(d), i) trajectories including two states, 
(2,1) and (3,2) (as indicated by thin arrows), and ii) trajectories including a state (2,1), 
but not including a state (3,2) (as indicated by broken arrows). However, the state (2, 1) 
is no longer the attractor although (3, 2) is still an attractor. Then, the diagram includes 
trajectories involving an attractor (3, 2) and those involving no attractors. In the former 
trajectories, the escape rate from (3, 2) is still independent of F. Therefore, u increases 
slightly with F. 

In the case where / < F(< I + V), the escape rate from the attractor (3,2) becomes 
one increasing exponentially with F, which is given through the transition (3, 2) (3, 3) 
as in Fig. 4(e). Then, the mean velocity of the driven particle along trajectories including 
the attractor (3, 2) (as indicated by thin arrows) increases exponentially with F, although 
that along any trajectories without the attractor (3, 2) (as indicated by broken arrows) are 
unaffected by F. Therefore, within this range of F, u increases steeply with F. 

Thus, for V < 1/3, u increases slowly steeply — > slowly — > steeply with the increase 
in F. This result is consistent to the F — u relations in Fig. 2(a). 

In cases with //3 < V < 1/2 or//2 < V , the manner of the F dependent changes of 
the transition diagrams proceed differently from the previous case. (We do not show the 
diagrams for such V in order to avoid the complications.) We can also explain the mechanism 
of the phenomena observed in each range of V via the same analytic method. 
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SUMMARY AND DISCUSSIONS 



In this paper, we investigated the transport phenomena of a simple nonequihbrium lattice 
gas system. We measure the mean velocity of the driven particle {u) and the coefficient of 
effective drag of this particle {r]). Under the certain ranges of parameters, we obtained the 
u - rj relations qualitatively similar to the relations which observed between the shear rate 
and shear viscosity coefficient in the shear-thickening polymer solutions jl^. 15,^^ or that 
of the shear-thickening colloidal suspensions!^, |2^. The mechanisms of these phenomena 
are explained by the transition diagram. 

Now, we consider the similarities between our results and those obtained recent studies 
for the rheology of soft matter systems in a little more detail. The soft glassy systems 
like polymer solutions or colloidal suspensions under shear have been studied extensively by 
experiments, molecular dynamics simulations, the analysis of the transient network model 
or the random spin model, and studies of the mode coup ling theory or an analysis of the 



pair distribution function 



flDlilHHIilC, 



21 



Mm, 



Among results of some experimental and theoretical studies for shear-thickening polymer 
solutions, it has been observed that the shear viscosity varies in the form, increase — >■ fast 
decrease — >■ slow decrease fast decrease, with the increase in the shear rateQ,Q|. Such 
a profile was obtained in the u — rj relation in our simulation under a certain condition (for 
the case I-ii in the phase diagram Fig. 3(d)). On the other hand, among results of some 
experimental and theoretical studies for shear-thickening colloidal suspensions, it has been 
observed that the shear viscosity varies in the form, decrease — increase decrease, with 



the increase in the shear rate 



21 



25 



Such a profile was also obtained in the 



u — 7] relation in our simulation under the other certain conditions (for the case Il-ii or III 
in the phase diagram Fig. 3(d)). 

Moreover, among results of recent studies, the power law [shear viscosity coefficient] ~ 
[shear rate]"'' has been observed in the shear-thinning regime. Herej5 ~ 2/3 has been 



given in some experimental studies, simulations and theoretical studies 



2^,|2^, 



24 



2611 while 



5 that varies depending on several conditions between and 1 has been obtained in the 



other experimental and theoretical studies 



21 



2^ . On the other hand, we obtained the 



following u — 1] relations in the rj decreasing regime for lower T; rj u 



0.5 



for smaller u 



and 7] u for larger u. 
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Thus, we expect our results to provide important hints to uncover the possible mechanism 
for several rheological characteristics of several soft materials. However, the relations of some 
properties, for example the exponents of the transport coefficients, the transport properties 
in cases with lower temperature, etc. are still unclear between our results and those obtained 
in several recent studies. Then, studies of extended models, including more lanes or particles 
in the space with more sites or continuous space, and the comparisons between such toy 
systems and either real systems or more realistic models of the polymer or colloidal systems 
represent important future issues. 

The author thanks to M. Sano, and M. Otsuki for useful discussions. This research was 
supported in part by a Grant-in- Aid for JSPS Fellows (10039). 
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FIG. 1: Illustrations of effects of potential and external field in each lane. 




FIG. 2: Mean velocity u as a function of F for (a) V = 0.25 with T = 0.05 or T = 0.07, (b) V = 0.4 
with T = 0.05 or T = 0.08 and (c) V = 0.6 with T = 0.08 or T = 0.1. 
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FIG. 3: Coefficient of effective drag as a function of u for (a) V = 0.3 with T = 0.05 (Times) 
and V = 0.25 with T = 0.05 (Plus), T = 0.07 (Star), T = 0.2 (Empty square) or T = 0.32 (Full 
Square), (b) V = 0.35 with T = 0.06 (Times) and V = 0.4 with T = 0.05 (Plus), T = 0.08 
(Star), T = 0.2 (Empty square) or T = 0.32 (Full Square) and (c) V = 0.55 with T = 0.08 and 
V = 0.6 (Times) with T = 0.08 (Plus), T = 0.1 (Star), T = 0.2 (Empty square) or T = 0.32 (Full 
Square), (d) Phase diagram of the system as functions of V and T. Dots and broken lines give the 
boundaries between states 
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FIG. 4: Typical transition diagram (left), and illustrations of attractors and rough estimated 

escape rates from them (right) for (a)F < V, (h)V < F < I - 2V, (c)I -2V<F<I-V, 
{d)I -V <F <I and (e)/ <F <I + V under V < 1/3. 
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